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1 . INTRODUCTION 


One method of improving upon the performance of rotating machinery 
components such as motors, generators, turbines, and high speed bearings 
and drills is to provide internal cooling using the rotating, non- 
capillary heat pipe concept [1]\ 

2 

The rotating, non-capillary heat pipe is shown schematically in 
Figure 1. It consists of a sealed hollow shaft, having a slight internal 
taper from one end to the other, and containing a fixed amount of working 
fluid. When the shaft is rotated at high speed about its longitudinal 
axis, the working fluid collects as an annulus at the large end. The 
diameter may be stepped at this end to provide a larger liquid reservoir. 
Heat added to this end of the shaft (evaporator) evaporates the working 
fluid, generating vapor which then flows axially toward the other end. 

Heat removed from this end of the shaft (condenser) condenses the vapor. 
The centrifugal forces accelerate the liquid condensate back to the evap- 
orator to complete the cycle. 

1.1 OBJECTIVE 

The overall objective of this research program is to understand the 
operating characteristics and to predict the heat transfer performance 
capabilities of the rotating, non-capillary heat pipe. Preliminary work 
under this program has been reported earlier [2]. This report presents 
the final results obtained, with some recommendations toward further work. 


^ Numbers in brackets indicate references listed in the Bibliography. 
2 

This device should perhaps be called a rotating, two-phase, closed 
thermosyphon. 
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2. ANALYTICAL PROGRAM 

The analytical program initially focused upon the understanding of 
the operation of the rotating, non-capillary heat pipe, including a study 
of the mechanisms of nucleate boiling and condensation heat transfer, and 
two-phase, annular, counter-current flow. The heat transfer limitations 
imposed upon the rotating, non-capillary heat pipe by (a) the critical 
nucleate boiling heat flux ("burnout"), (b) the vapor-condensate entrain- 
ment velocity ("flooding"), (c) the sonic vapor velocity, and (d) the film 
condensation thermal resistance were estimated and compared. It was con- 
cluded that the performance of these devices would be controlled primarily 
by the overall resistance to heat transfer in the condenser, including 
the Internal film condensation resistance, the conduction resistance across 
the condenser wall, and the outside surface heat transfer resistance 
[2, 3, 4]. The rotating non-capillary heat pipe was then modelled as a 
rotating, truncated cone, and a Nusselt-type film condensation analysis 
was performed. (Figure 2) 

2.1 FILM CONDENSATION MODEL 

Laminar film condensation with rotation has been analyzed for a vari- 
ety of geometrical conditions. In 1959, Sparrow and Gregg [5] obtained a 
similarity solution for film condensation on a rotating disk situated in 
a quiescent vapor. They formulated the problem as an exact solution of 
the Navier-Stokes and energy equations, and arrived at numerical results 
for a wide variety of test fluids. Later, in 1961, Sparrow and Hartnett 
[6] showed that the disk solution could be modified to apply to film 
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condensation on the outside of rotating cones which are not too slender. 

In 1963, Singer and Preckshot [7] studied condensation on the outside of 
a rotating shaft. They reported the presence of three condensate flow 
regimes, which depended on shaft speed, and theoretically estimated the 
heat transfer during two of the regimes. Chato [8], in 1965, studied 
film condensation in a variable gravity field during rotation of a 
condensing thermosyphon. He included the effect of surface shear between 
the liquid condensate and the moving vapor. 

During this work. Ball back [3] performed a Nusselt-type analysis 
for film condensation on the Inside of a rotating, truncated cone of 
half angle 4 as modelled in Figure 2. Ball back made all the familiar 
assumptions used in Nusselt's theory of laminar film condensation, includ- 
ing negligible subcooling and momentum changes in the condensate film, 
no interfacial shear between the condensate and vapor, and pure conduc- 
tion within the condensate film. In addition, he assumed that the cen- 
trifugal acceleration is much larger than the normal acceleration of 
earth gravity, the thickness of the film is much less than the radius of 
the condenser wall, and the vapor space is essentially stagnant. In 
order to analytically solve for the condensate film thickness, however, 
Ballback further assumed that the slope of the condensate film, ds/dx, 
is much less than tan $. This last assumption is not valid for very 
small half-cone angles $. He analytically solved for the condensate film 
thickness and heat transfer coefficient, and his results were later confirmed 
as a special case of the work of Dhir and Lienhard [9] who studied laminar 
film condensation on axisynmetric bodies in non-uniform gravity fields. 

Daley [4] reasoned that most practical applications of the rotating, 
non-capillary heat pipe would involve very small half-cone angles. He 
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therefore modified Ballback's work by removing the restriction that 
da/dx is much less than tan <*>. He also improved upon the analysis by 
including the thermal resistances due to the condenser wall and to the 
outside surface cooling mechanism. Using the velocity profile for the 
condensate developed by Ball back, Daley arrived at a non-linear second 
order differential equation for the film thickness, 6 [2]. He numeric- 
ally integrated this expression using a Runge-Kutta-Gill numerical 
integration scheme. To start the integration, he used the conditions 

that 6 = 6. and d<s = tan $ at x « 0. It was postulated that 6, was a 
1 dx 1 

function of the over- fall condition at the step located at the condenser 

exit. Daley arrived at an expression for the minimum film thickness, 

6 min a * the condenser using the approximation that the flow over 

the step approximated the free fall condition for open channel flow. 

Having found the film thickness from the second order differential equa- 
tion at x * L, the condenser exit, Daley compared this value with a^at L. 
Using this technique, Daley was able to find the heat flow out of the con- 
denser as a function of RPM for 0° half cone angle only. As the cone 
angle was increased to 1°, however, the numerical integration of the 
highly non-linear second order differential equation became extremely 
sensitive to the initial value of the film thickness. Consequently, no 
solutions were found for half cone angles greater than zero. 

The Nusselt type film condensation analysis performed by Ballback 
and Daley was extended by Newton [10] to include the effects of vapor 
pressure drop and Interfacial shear. Equations for momentum, continuity 
and energy were developed for an infinitesimal fluid element in the con- 
densate, and, assuming one dimensional turbulent flow, momentum and 
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continuity equations were developed for Isothermal vapor flow. The details 
of his analysis are now Included for completeness. 

2.2 ANALYSIS OF NEWTON ClOl 
Condensate Momentum Equations 

Assuming that the momentum changes within the condensate film are 

small, a static force balance on a differential liquid element In the 

x-dlrectlon (See Figure 3) yields: 

0 • h- - ifi. + p, w 2 r sin * n 1 

ay ax ' 

where 

t * local shear stress 
p » local pressure 

x • co-ordinate measuring distance along 
the condensing surface 

y » co-ordinate measuring distance normal 
to the condensing surface 

r * local radius in the condensate 

4 3 half-cone angle 

Pf * density of the condensate film 

(i> * angular velocity 

Similarly, a force balance in the y-directlon yields: 

2 

o * - aj> - p u r cos$ (2) 

a / T 

In the above equations, the local radius at any point in the liquid con- 
densate may be expressed as: 

r(x,y) * Ro ♦ x sin$ - y cos$ (3) 

where Ro 3 minimum wall radius in the condenser section. 
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Newton [10] substituted Equation (3) for r(x,y) into (2) and integrated 
from p to p y and y to 4 (the local condensate thickness) to obtain the 

local pressure at any point within the condensate: 

2 

p(x,y) = p v (x) + p f to cos* {(Ro + x si n4>) (6 - y) 

- 1_ cos* (6 2 -y 2 ) } (4) 

2 


Differentiating Equation (4) with respect to x, substituting this result 


into Equation (1), and integrating from y to 6 and t to -t y gives: 


t = 


dp v (6 -y) - t + p u 2 [Ro + jsi n«|> - 5cos*][sin* - cos* d6]( 5-y) 
W v f dx 


a iLM. (xjJTl 

sly 


(5) 


where 

u f = viscosity of the condensate film 
u(x,y) = velocity of the condensate film 
Another integration from 0 to y and from 0 to u gives an expression for 
the velocity: 

u(x,y) 3 - J_ 5v (y« - xl) - JV y 

Uf dx 2 

2 

+ p f“ [Ro + x sin* - 5cos*] [sin# - cos* d6] (y6 - y^) (6) 

dx 2 


Continuity Equation 

Conservation of mass, with condensation occurring, requires that 

dWj. _ dw v 
dx dx 


(7) 


or 



(7a) 
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where 


w- 3 mass rate of flow of condensate 

« -A 


a jp f u(x,y) 2irr(x,y)dy 

(8) 

w y 3 mass rate of flow of vapor 


3 p y trR 2 (x) v (x) 

(9) 

and 


R(x) 3 Ro + x sin<|> - 6 cos$ 

(10) 

Energy Equation 



An energy balance on the differential element may be written 
as: 


2tt(Ro •>• x sinji) (T t - T°°) (11) 

6 . t . 1 

1 

where 

T $ = saturation temperature of the vapor 
T« * ambient coolant temperature 
t 3 thickness of the condenser wall 
6 = film thickness of the condensate 

k^ = thermal conductivity of condensate film 
k 3 thermal conductivity of condenser wall 

W 

h 3 external heat transfer coefficient 
hfg 3 latent heat of vaporization 
q = heat transfer rate 

Note that the above differential equation is only valid for small half- 
cone angles (where the condenser geometry is not too different from cyl- 
indrical) and for thin condenser walls. 
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Vapor Momentum Equation 

Using the differential element shown in Figure 4, and neglecting 
the momentum contribution due to the condensing vapor, a balance of 
forces may be taken in the x-direction and set equal to vapor momentum 
changes: 

-A, ^ + V, " 02) 

where 

P v (x) = pressure of the vapor 
v(x) * velocity of the vapor 
A y (x) » cross sectional area of vapor element 
* ttR(x ) 2 

S y (x) = perimeter of vapor element 
= 2 -jtR(x) 

Expanding the right hand side of Equation (12), and using Equation (9), 
gives the following equation for the vapor pressure gradient: 

-2p v dv - 2v 2 p dR + 2 _Tv (13) 

dx dx R dx R 


The first two terms in Equation (13) represent the momentum effects of 
the vapor while the last term represents the usual wall -friction effect. 
The interfacial shear stress may be written 



f p v v(x ) 2 


(14) 


For low condensate velocities, according to Bergelin, et al [11], the 
two-phase friction factor may be approximated by the friction factor used 
in smooth tubes. For simplicity, the assumption was therefore made 
that the effect of the liquid film on the condenser surface is negligible. 


8 


and that 


where 


f = constant 

R .' 


( 15 ) 


R g = vapor Reynolds number 
a 2R(x)pyV(x) 
y v 

The constant and the exponent in Equation (15) would, of course, depend 
on the Reynolds number range, and are readily available in heat transfer 
textbooks [13]. 

A more refined analysis could be performed, using the Lockhart-Martinell i 
correlation [12] for two-phase frictional pressure drop, but since the 
friction effects were considered to be small in this study, the refined 
analysis was not warranted. 

The above system of equations can be reduced to two first-order non- 
linear differential equations with unknowns v(x) and 5(x). 


dv 

37 


and 


2(Ro + x sin<j>) (T $ - T®) 




sin$ - cos 4 > d6 
37 


(16) 


dx 

where 


P] (x) ^ P 2 

(x)' 

[ “f R < x > 2 d»l 

p v f v(x) 2 R(x) 

T~ . 

1 PvV(x) L 2p f 

- 2Pi(x)37J 

Pi(x) cos$ 

p f u 2 R(x) 2p v(x) 2 



L 

+ KTTj 


= (Ro + x si n<|») 6 (x) 3 

- ^ COS4> s(x) 4 



? 2 (x) 3 (Ro + x sin<i>) . 6(x) 2 - cos$ 5(x) 3 

~cL~ -r- 
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A Runge Kutta numerical integration was used to solve these equations 
with an_IBM 360 digital computer. It was hoped that-by using two-first- 
order differential equations rather than a single second order equation 
as Daley did, that the sensitivity to the initial value of 6 would be 
reduced. The integration was started by assuming initial conditions, 

6 = 6. and v = 0 at x = 0. With these initial conditions = tan<j> at x=0. 

1 dx 

Figure 5 is a schematic drawing representing the calculated 6 (x) for dif- 
ferent initial values 6^. Note that if 6. is too small the calculated 6 
becomes negative before the condenser overfall (physically impossible solu- 
tion). If 6^ is too large the calculated 6 grows at a constant slope equal 
to tanc(> (physically unrealistic). With the proper choice of 6^ the calculated 
6 reaches a maximum along the length of the condenser and then tapers off to 
a value close to zero at the overfall. The exact value of 6 at the overfall 
is not known. It was observed however that the calculated value of 6 at the 
overfall was approximately 0.1 of 6^ , and the heat transfer rate was insensi- 
tive to an exact value at the overfall . This result agrees with Nimmo and 
Leppert [14] who experimentally found for laminar film condensation on a 
finite horizontal surface (with proper overfall drainage) that the heat trans- 
fer rates were practically independent of the end conditions. 

Newton obtained results for half-cone angles of 0°, 0.1°, 0.2° and 0.3°. 
As the half-cone angle increased, the solution of the differential equations 
became extremely sensitive to the initial values of 6^. He also noted that 
his results for 6(x) snowed that (except near the ends of the condenser) 

<< tan<j>. Newton, therefore, made this approximation, and simplified the 
model further using an order of magnitude analysis of the terms in Eqs. (16) 
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and (17). With these simplifications, Eqs. (16) and (17) were reduced to: 

^ - 2(Ro * xsin») (Tt - T»,) - 2v(>)sint (18) 

p»|«M ♦i*lK*(x ) 2 R(x) 

kf *w ' 


6(x) 3 R(x)sin * (p f w 2 R(x) + 2 p y v(x) 2 ) - 6(x) 2 P y fv(x) 2 R(x) 


R(xT~ 


- p„v(x) 


Equation (18) was numerically integrated using the initial condition 
that 6^* 0 and v * 0 at x = 0, and using the cubic equation for 6 (x) * 
Equation (19). After the integration was completed for 6(x) and v(x), 
then the energy equation was integrated over the length of the condenser 
to find the total heat transfer rate of the heat pipe. Newton's results 
for both Equations (16) and (17) and also for Equations (18) and (19) are 
tabulated in Table I. 


Table I 

Total Condenser Heat Flow Rate q (BTU/hr) 


* 

Sol'n. to 

Equations (16) & (17) 

Sol'n. to 

Equations (18) A (19) 

0° 

8,505 


0.1° 

10,082 

9,836 

o 

CM 

• 

O 

11,173 

11,033 

0.3® 

11,904 

11,769 


Notice that the approximate solution is slightly conservative, and becomes 
more accurate as the half-cone angle increases. Since most practical 
designs will involve half-cone angles at least as large as 0.5° to 1.0°, 
the approximate solution can, therefore, be used with little error. 


11 



2.3 HEAT PIPE PARAMETRIC STUDY 


A parametric study of the heat pipe operating characteristics was 
performed using this approximate solution [15]. In an attempt to estab- 
lish these parametric relationships, a preliminary dimensional analysis 
was performed on the system of approximate equations which led to the 
functional relation 

9 s ^/kJLLM, “Bs. 2 , Kw, Jl ,1152., , — Hit — ^(20) 

k^L(T s ~T oo ) \ Ro k^ Ro k^. c^AT/h^g ^f^o^Vg 

where q is a non-dimensional form of the heat transfer rate 

k^-T.) 

of the heat pipe. The independent operational variables were separated 
into broad categories of (1) working fluid properties, (2) operating con- 
ditions, and (3) heat pipe geometry. From these categories convenient 
parameters were selected to be studied. 

Woodard [15] obtained results for a stainless steel heat pipe, with 
a 1/16 inch thick condenser wall, 9 inches long. Working fluids were 
water, ethyl alcohol and freon 113; condenser half-cone angles were 0°, 

1°, 2° and 3°; condenser minimum internal radii were 0.125", 0.25", 

0,5", 0.75", 1.0" and 2.0"; RPM was varied in increments of 600 RPM 
from 600 to 3600 RPM and values for the exterior heat transfer coefficient 
were 50, 100, 500, 1000, 3000 and ». The complete results are contained 
in Woodard's thesis [15]. 

Effect of Working Fluid 

Figure 6 shows the heat transfer response of the rotating heat pipe 
to the overall physical properties of the working fluid in use. Water, 
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in the range of operation studied, is obviously the most satisfactory fluid. 
Its high latent heat of vaporization and high thermal conductivity relative 
to that of ethyl alcohol and freon 113, would explain this difference. For 
each working fluid, the heat transfer performance is improved with heat pipe 
operating temperature (and pressure). 

Effect of Rotational Speed 

Ball back's analysis indicated that the heat transfer capability of 
the heat pipe is proportional to the square root of the RPM. This behavior 
is weakened in Newton's analysis by the additional thermal resistance due 
to the condenser wall and the external heat transfer coefficient. However, 
analysis of the trend presented in Figure 7 indicates that the heat pipe 
response approaches this theoretical relationship, q * Ju~ as the ex- 
terior heat transfer coefficient becomes very, very large. 

Effect of Exterior Heat Transfer Coefficient 

Figures 7 & 8 suggest the dependence of the heat pipe capability 
on the condenser cooling mechanism. It is easily seen that the exterior 
heat transfer coefficient dominates the heat pipe capability when the 
exterior heat transfer mechanism is poor. 

Effect of Heat Pipe Geometry 

Figure 9 shows that the heat transfer rate increases rapidly with 
increasing half cone angle, from o° to 1°, for a nominal RPM and heat 
transfer coefficient. This behavior then approaches an approximately 
linear response to half cone angles of 3 degrees. Figure 9 also shows 
that there is a significant heat transfer rate even when there is no in- 
ternal taper to the condenser (when <f> = o°). 
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The change of heat transfer capability with change in internal radius 


is i llustra ted by Figure 10. The behavl or_of_.the_hea t_-transf er_ra te i s 

related to the radius by two mechanisms, the change in condenser area and 
the change in centrifugal acceleration at the condensate film. 

2.4 SUMMARY OF ANALYTICAL RESULTS 

1. The approximate numerical solution [using Equations (18) and 
(19)] is valid for most practical heat pipe geometries. 

2. In the temperature range from 80° F to 240° F, water appears to 
be the best working fluid. 

3. The heat transfer capability of the heat pipe is dominated by 
the exterior heat transfer mechanism when this mechanism leads to low 
exterior heat transfer coefficients. 

4. There is a significant condenser heat transfer rate even with a 
0° half cone angle heat pipe. After an initial rapid rise with half 
cone angles from 0° to 1°, the heat transfer rate continues to increase 
approximately linearly with increasing half cone angles. 

5. The heat transfer rate of the heat pipe approaches a function 
of the square root of RPM for efficient exterior heat removal mechanisms. 

6. A change in the minimum internal condenser radius, Ro has a 
large effect on heat pipe heat transfer. 
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3. EXPERIMENTAL PROGRAM 


3.1 DESCRIPTION OF EQUIPMENT 

The main components of the rotating, non-capillary heat pipe are 
grouped into the evaporator, condenser, auxiliary equipment, and instru- 
mentation. Figure 11 is a schematic diagram of the test apparatus. A 
cross-sectional drawing of the assembled heat pipe is pictured in Figure 
12, and Figure 13 is a photograph of the experimental facility. 

Evaporator 

The evaporator is a 3.125-inch inside diameter stainless steel 
cylinder, 5.90-inches long. One end is flanged to an outside diameter 
of 5.906-inches to accommodate the condenser and to support a large, 
single row precision ball bearing. The other end is flanged to an out- 
side diameter of 4.50-inches to accommodate two pyrex glass end windows. 
The inner glass window is fused to a i /4-inch diameter glass filling tube 
which passes through a hole in the outer window. The inner window presses 
against a neoprene o-ring and is separated from the outer window by a 
compressed fiber gasket. Both windows are held in place by a stainless 
steel end cap. (See Fig. 12.) 

The evaporator is helically wound with an 11-gage Chromel-A heater 
wire in a 3/16-inch outside diameter Inconel sheath. The heater coils 
are bonded in place using a gold-nickel brazing alloy in a controlled 
hydrogen atmosphere. They are coated with a thin layer of Sauereisen 
cement and are packed with asbestos insulation to reduce radial heat 
losses. In addition, a 1/16-inch wide radial groove is machined into 
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the evaporator wall on either end of the heater element to reduce axial 

heat Josses. Electrical power. to_ the heater _is_ passed -by .a. graphite 

brush assembly through bronze collector rings. Power is furnished to 
the heater coils by a 120 volt DC motor generator set with an automati- 
cally regulated voltage output, and is controlled using five resistive 
load banks in series with the heater coils. 

Condenser 

The condenser is a 10.0-inch long stainless steel, truncated cone 
with a 1.0 degree half-cone angle and with a 1/1 6-inch wall thickness. 

The large end of the condenser is flanged to bolt to the evaporator and 
has an inside diameter of 1.81-inches. The small end is machined into 
a cylinder 1.46-inches inside diameter, by 1.50-inches long. It is 
flanged to accommodate a cylindrical end plug. Each flanged joint is 
sealed with a teflon-coated, metallic o-ring. 

The stainless steel condenser end plug is hollowed out to allow the 
pressure transducer arm to be passed through its inside face, and to 
allow the drive shaft to be threaded and keyed in place. The outer end 
of the plug is machined down to four flat sides, each 7/8-inch wide and 
7/8-inch long to accorranodate phenolic thermocouple junction boards. 
Auxiliary Equipment 

The auxiliary equipment are grouped into the drive assembly, test 
stand, spray cooling assembly, safety shields, and filling system. 

Drive Assembly 

The drive assembly consists of the drive shaft, support bearing, 
pulley and variable drive motor. The drive shaft is a 3/4-inch diameter 
stainless steel cylinder which is hollowed out to allow the instrumentation 
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leads to be connected to the slip-ring unit. The outside diameter of the 
shaft is stepped in several places to thread into the condenser end plug, 
to accommodate a 2.65-inch diameter drive pulley, and to support a double 
row, angular contact bearing. Each end of the shaft is internally threaded, 
The pressure transducer arm is screwed into one end and the slip-ring coup- 
ling into the other. The shaft is screwed securely into the condenser end 
plug and is keyed in place. Torque is applied to the shaft by a V-belt 
using a 2HP, 3 phase variable speed motor. The motor is capable of speeds 
from 450 to 4500 RPM and is equipped with a magnetic disc brake and an 
electric remote control unit. 

Test Stand 

The test stand is designed so that the heat pipe can be rigidly sup- 
ported with its longitudinal axis ranging in position from 0 to 90 degrees 
from the horizontal. Both the heat pipe and variable drive motor are 
mounted on steel bed-plates connected by a 4-inch diameter iron pipe. The 
iron pipe is held in place by three 2-inch thick clamps which are supported 
two feet off the ground by a 1/4 inch steel plate. As the heat pipe is 
rotated in any position, the drive motor rotates with it and the entire 
assembly is then securely held in place. 

S pray Cooling Assembly 

Cooling of the condenser section is done by spraying tap water on 
the outside of the condenser during rotation. The condenser section is 
enclosed by a 9-inch long 13-inch diameter well- insulated stainless steel 
cylinder. The ends of this cylinder have Garlock rubber seals separating 
them from the rotating condenser section. The condenser cover is slid into 
place over the condenser section through its open ends and is then bolted 
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in place to the support plate. The cover has an 8-inch x 6-inch door 
on it to make adjustments-or repairs on the heat _pipe without removing 
the entire cover. The cover has eight spray nozzles. Each nozzle is 
mounted at the same axial position but is spaced 45° apart on the circum- 
ference. Various nozzles may be used to give a range of water droplet 
sizes from 300 to 600 microns. The cooling water flows from copper feed 
lines to a mixing tube soldered onto the top of the condenser cover. The 
coolant flows from the mixing tube through insulated plastic tubing to 
each of the spray nozzles. It drains from the bottom of the condenser 
cover to another mixing tube before being dumped into the building drain 
lines. 

Safety Shields 

To insure safe operation, the entire heat pipe assembly is surrounded 
by 1/8 inch thick stainless steel shielding which easily bolts to the 
support plate. 

Filling System 

The filling system consists of a glass filling manifold, mercury 
manometer and vacuum system. The glass manifold and manometer are shown 
schematically in Figure 14. A mechanical vacuum pump and air-cooled 
oil diffusion pump are used to evacuate the system. 

Instrumentation 
Temperature Measurements 

Eight Kapton insulated copper-constantan thermocouples are used on 
the rotating heat pipe. Four thermocouples are manufactured with a 1/16 
inch diameter stainless steel sheath. Each of these is placed within a 
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1/16 inch diameter well drilled 3 1/2 inches deep into the evaporator wall 
at different radial positions to monitor the radial heat transfer into the 
evaporator. One thermocouple is located 0.066 inch radially outward from 
the inside of the evaporator surface; another is located 0.188 inch radi- 
ally outward; the remaining two are each located 0.125 inch radially out- 
ward, but are spaced circumferentially 180° apart to monitor temperature 
uniformity. Each radial location has an estimated uncertainty of ±0.005 
inches. Two thermocouples, sealed in stainless steel sheaths, are suspen- 
ded in the evaporator vapor space to measure the saturation temperature of 
the vapor. Two additional thermocouples are spot welded to the outside of 
the condenser at different axial positions to monitor the heat transfer 
through the condenser. The leads from all these thermocouples are brought 
along the sides of the condenser and through the small rear condenser 
flange to the phenolic junction boards on the condenser end plug. From 
the junction boards, the leads go inside the shaft to a mercury slip-ring 
unit, rated at less than 10 micro-volts noise at speeds to 4000 PPM. The 
slip-ring output is fed to a Hewlett-Packard 201 OC Data Acquisition Sys- 
tem which has an accuracy of ±0.5 microvolts. 

All evaporator wall thermocouples were calibrated with an uncertainty 
of ±0.5°F using the triple point of water (32.0°F), the boiling point of 
water (212. 0°F at 14.696 psia), the melting point of tin (449. 4°F), and 
the melting point of lead (621. 3°F). The other thermocouples were cali- 
brated with an uncertainty of ±0.2°F using the triple point and the boil- 
ing point of water. 

The temperature of the condenser cooling water is measured before and 
after spraying the test section by quartz crystal thermometers mounted in 
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well -insulated mixing tubes located above and below the condenser spray 
cooling cylinder. These thermometers have an accuracy of better than 
±0.1 °F. 

Pressure Measurements 

The saturation pressure of the vapor is measured by a semi-conductor 
pressure transducer mounted in the vapor space at the small end of the 
condenser. The transducer is threaded on a 1 /4-inch diameter, 2-inch 
long arm which extends from the shaft through the end plug and into the 
vapor space. The leads of the probe pass through the transducer arm, 
through the shaft to the junction boards, and out to the slip-ring unit. 

The pressure transducer was calibrated against a Wal lace-Tiernan 
standard pressure gage (with an accuracy of 0.03 ± 0.01 psia) for pressures- 
greater than atmospheric, and against a mercury manometer (with an uncer- 
tainty of ±0.02 psia) for pressures less than atmospheric. The calibra- 
tion was performed at different operating temperatures by placing the 
transducer within a controlled temperature oven with a quartz crystal 
thermometer. 

Coolant Flow Rate Measurements 

The coolant flow rate is measured by passing the coolant through a 
calibrated flow rotameter (with an accuracy of ±2%) before it enters 
the spray nozzles. 

Rotational Speed Measurements 

The rotational speed of the heat pipe is measured using a geared 
fly-wheel and a variable reluctance transducer whose output frequency is 
recorded using a Systron counter. An electronic strobe light is used 
as a check. 
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Electrical Power Measurements 

Electrical power to the heater coils is measured by an ammeter and 
voltmeter. 

3.2 EXPERIMENTAL PROCEDURES 

The experimental procedures consist of cleaning and filling pro- 
cedures and operational procedures. 

Cleaning and Filling Procedures 

In order to obtain film condensation in the heat pipe condenser when 
water was the working fluid, the interior surface was prepared to insure 
that proper wetting occurred. In addition to wetting, it was desirable 
to remove as much of the non-condensible gases which were present, and 
to fill the heat pipe with the same amount of fluid each time. The above 
goals were accomplished by performing a detailed cleaning and filling pro 
cedure as described below. 

a. The heat pipe was placed in the vertical position and mounted 
to the fill equipment. 

b. With end windows and "0" rings removed, the heat pipe was filled 
with a chlorinated hydrocarbon degreasing compound such as Tri- 
chloroethylene. 

c. The pipe was then drained and filled with distilled water. 

d. Steps b. and c. were repeated. 

e. The pipe was then flushed three times with alcohol. 

f. Step c. was repeated. 

g. The pipe was filled with a warm mixture of 8 parts distilled 
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water, 2 parts ethyl alcohol, 2 parts 50% NaOH solution and 1/2 part 30% 

solution. 

h. Itwas then thoroughly rinsed wi th distilled water, while noting- 
any signs of non-wetting. 

i. The pipe was kept filled with water until ready for filling. 

When the heat pipe was ready to be filled, it was prepared in the follow- 
ing manner. 

a. All the water was drained from the pipe, and the end windows 
were mounted in place. 

b. The glass filling tube was connected to the filling manifold, 
and a vacuum was pulled on the test section for about 1/2 hour. 

c. Approximately 150 ml. of de-gassed, distilled water was then 
slowly bled into the heat pipe. 

d. The glass fill tube was then sealed under vacuum. 

Operational Procedure 

Having completed the filling procedure, the heat pipe was ready to 
run at a pre-determined RPM. About 10 drops of lubricating oil were added 
to the oil fill tube at the top of the large bearing. The cooling water 
lines were turned on so that the rotameter was at approximately 30% flow. 
The drive motor was turned on and the heat pipe was slowly brought up to 
the desired RPM. Once a liquid annulus was formed in the evaporator and 
the desired RPM was reached, the DC motor generator was turned on and 
power was added to the evaporator. After thermal equilibrium was reached 
for a particular power setting, data was taken by the instrumentation de- 
scribed earlier. Shut-down was accomplished by first shutting off the 
current to the heater, and then decreasing the RPM by use of the remote 
control unit. 
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3.3 REDUCTION OF DATA 

The condenser heat transfer rate was experimentally determined using 
the measured mass rate of flow and the overall temperature increase of 
the cooling water in Equation (21). 

q = wc p AT (21) 

A correction was made to this result, however, to take into account 
frictional heating effects of the bearings and seals, and viscous dissi- 
pation effects within the cooling water during rotation. To measure 
these effects the heat pipe was rotated at various rotational speeds with 
cooling water flowing through the spray nozzles, but with no input power 
applied to the evaporator. In each case, the mass rate of flow of the 
cooling water was measured along with the overall temperature increase 
of the cooling water across the condenser box. These measured values 
were used in Equation (21) to calculate the heat removed due to friction- 
al effects at various rotational speeds. This calculated frictional heat 
transfer rate was subsequently subtracted from the total condenser heat 
transfer rate for each selected rotational speed. 


23 



4. EXPERIMENTAL RESULTS 


4.1 EVAPORATOR PERFORMANCE 

During the initial operating runs, the evaporator performance was 
studied using visual observations and the measured evaporator wall temper- 
ature profile. It was noted, for instance, that at high rotational speeds, 
the number of active nucleation sites decreased and bubble departure dia- 
meters were smaller. This result agrees with the observations of Gray & 
Marto [16, 17] during nucleate boiling within a cylindrical, rotating 
boiler. In addition, at low saturation pressures (corresponding to low 
input powers), the nucleate boiling action was violent, with large vapor 
bubbles and a frothy liquid-vapor interface. On the other hand, at high 
saturation pressures, the boiling action was more controlled, with smaller 
bubbles. 

Temperature and pressure data were recorded, and the measured temper- 
ature within the vapor space compared to within 1°F of the saturation tem- 
perature calculated from the measured vapor pressure. Moisture within the 
pressure transducer, however, prevented its further use due to unpredict- 
able shifts in its zero reading. 

Using the measured radial wall temperature profile in the evaporator, 

and assuming one-dimensional conduction, the heat flux and the inside wall 

surface temperature were calculated at the axial midpoint of the evaporator. 

Figure 15 presents the calculated heat flux versus (T - T ) for a 700 RPM 

w s 

run (corresponding to 22 G's at the test surface). Notice that as the 
saturation temperature in the vapor space T g increases, due to an increase 
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in the operating pressure with applied power input, the wall superheat 
decreases. The heat transfer coefficient therefore increases with applied 
heat load and operating pressure. This result is in agreement with the 
well-known pressure effect upon saturated nucleate pool boiling [13]. The 
data of Ponter & Haigh [18] at 1G, which shows this pressure effect, is 
plotted in Figure 15. In addition, the data of Gray, Marto & Joslyn 
[16], at 25 G's is shown for comparison. The first six experimental 
data points are in reasonable agreement with the other investigations. 

Data point 7 however, appears to be in error, perhaps due to the large 
uncertainty in the extrapolated inside wall surface temperature at high 
power inputs. 

4.2 CONDENSER PERFORMANCE 

A great deal of difficulty was initially experienced in trying to 
obtain good wetting of the condensate upon the stainless steel condenser 
surface [10]. The poor wettability of the condensate led to break-up of 
the condensate film with rivulet formation. After attempting several 
unsuccessful cleaning procedures, the one described above was tried and 
good wetting was obtained. The condensate then appeared as a mirror smooth 
film on the inside of the condenser surface. No ripples, rivulets or waves 
were observed. 

Effect of Dropwise Condensation 

Figure 16 compares the condenser heat transfer results at two rotational 
speeds of 700 and 2100 RPM. The heat transfer rate was calculated from the 
measured mass rate of flow of the cooling water, together with the measured 
temperature increase across the coolant mixing box, as noted earlier. Data 
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was taken during smooth film condensation (with good wetting) and also with 
dropwise condensation (with poor wetting). The dropwise condensation mode 
was promoted by wiping the inside of the condenser surface with an oily 
cloth. From this data, it is evident that the performance of the condenser 
can be improved with higher RPM, and with dropwise condensation. Note, 
however, that the influence of rotational speed is felt more strongly in 
the film condensation mode than in the dropwise condensation mode. The 
effect of high centrifugal acceleration may tend to flatten out the conden- 
sate drops, rendering them less effective in their mechanism of heat re- 
moval. This behavior (i.e. the influence of high centrifugal acceleration 
upon dropwise condensation) should therefore receive further study. Also, 
at a given RPM, it appears that the increase in the heat transfer capability 
of the heat pipe in going from film condensation to dropwise condensation 
is more pronounced at lower saturation temperatures, corresponding to lower 
operating pressures. This behavior also warrants additional study. 

Effect of Non-Condensible Gases 

In Figure 17, the heat pipe condenser performance is compared for two 
different filling procedures. For Run 1, the heat pipe was evacuated for 
10 minutes to 19 mm Hg before filling and sealing. On the other hand, for 
Run 3, the heat pipe was evacuated for 30 minutes to 4.5 mm Hg prior to fill- 
ing and sealing. The difference in behavior for these two runs is therefore 
due to varying amounts of non-condensible gas trapped within the heat pipe. 
Film condensation occurred during both runs. Notice that at 700 RPM, as 
well as 2100 RPM, better performance occurred during Run 3 which presumably 
had a smaller amount of non-condensible gas present. This result agrees 
with the well-known influence of non-condensible gases on condensation heat 
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transfer [13]. In general, the non-condensible gases blanket the condenser 
surface and prevent the migration of the vapor molecules toward this cold 
surface, reducing the heat transfer. In a conventional, capillary-wick heat 
pipe (with no rotation) the vapor movement pushes all the non-condensible 
gases down to the far end of the condenser where a vapor-gas interface may 
exist. This interface is not sharp however because of diffusion effects of 
the gas molecules through the vapor [19]. In addition, gravitational ef- 
fects may also lengthen this interface. Within the rotating heat pipe, the 
centrifugal acceleration may force the heavier non-condensible gas (air) 
toward the evaporator, spreading out the vapor-gas interface along the 
entire condenser surface. Such an effect should be more important at low 
operating pressures where the steam density is very small compared to the 
density of air (see Figure 17). 

Condenser Surface Temperature Profiles 

During Runs 1, 2 and 3 an unexpected longitudinal temperature gradient 
was measured on the condenser surface using the thermocouples mentioned 
earlier. A marked increase in condenser surface temperature occurred 
toward the evaporator section of the heat pipe. Additional copper- 
cons tantan thermocouples were therefore spot welded to the condenser sur- 
face to provide for a better temperature profile. Runs 4, 5, 6, and 7 
were made while these more complete temperature profiles were measured. 

[20]. Figures 18 & 19 show these profiles for Run 4 at 700 and 2100 RPM 
and at different input powers. The exact cause of this non-uniform wall 
temperature is not known, however, this may be indicative of (a) non- 
condensible gases in the vapor space, (b) a poor coolant spray distribution 
upon the condenser surface (see Fig. 12), (c) a longitudinal heating effect 
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due to conduction from the electrical heater or from frictional heating 
in the large bearing and Garlock seal. 

4.3 COMPARISON OF THEORY AND EXPERIMENT 

Figures 20 through 23 compare the experimental heat transfer rates 
to the predicted theoretical rates. Some typical experimental error 
limits are included. 

In arriving at the theoretical results, several modifications were 
made to Eqs. (11 )&(1 8) . Since the condenser wall surface temperature was 
measured during these runs, it was assumed to be known in the theoretical 
model. The outside surface heat transfer coefficient was an unknown and 
was therefore not used. Equation (18) was consequently modified as fol- 
lows: 



Notice that the total thermal driving force becomes T s - T w (x) and the 
total resistance includes the condensation resistance and the wall conduc- 
tion resistance. The outside heat transfer resistance is unnecessary 

since T (x) is known. Equation (11) was modified in a similar way. 
w 

To arrive at an analytical expression for T w (x), the measured wall 
temperature profiles (See Figures 18 and 19) were used. A cubic poly- 
nominal approximation to the data was made using a least squares technique 
[20]. The resulting cubic polynominal was used for T w (x) in Equation (22), 
and together with Equation (19), a numerical integration was performed. 
During the integration, condensate properties were evaluated at a film 
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temperature equal to the arithmetic average of the saturation temperature 
and the local inside wall temperature. 

Figure 20 shows the comparison for Run 4 with water at 700 & 2800 
RPM. The theoretical results are approximately 20 percent higher than 
the experimental data at 700 RPM, and about 10 percent higher at 2800 
RPM. Such agreement is reasonable during film condensation heat transfer. 
Similar results are evident in Figure 21 for Run 5 with water. Figure 22 
compares the theory to the experimental data for ethyl alcohol, and Fig- 
ure 23 is for freon-113. Notice that for water, the theoretical results 
are higher than the experimental data, but for ethyl alcohol and freon, 
the theoretical results are less than the experimental data. This opposite 
trend was also recently found by Lee & Mital [21] during operation of a 
vertical, stationary closed thermosyphon. They found that their theory 
overestimates the heat transfer rate for water and underestimates the 
rate for freon-11. The reason for this discrepancy is not clear at this 
time. It may be due to differences in vapor densities between water, 
alcohol and freon 113, leading to differences in non-condensible gas 
migration, or to differences in vapor velocity effects within the theoret- 
ical model. In general, however, the agreement is reasonable and the 
approximate theory can be used to predict rotating heat pipe performance 
to within +20%. 
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5. CONCLUSIONS AND RECOMMENDATIONS 


5.1 CONCLUSIONS 

Based upon the experimental results obtained for the range of vari- 
ables covered during this investigation, the following conclusions can 
be made: 

1. The rotating, non-capillary heat pipe can be used effectively to 
transfer large quantities of heat in rotating systems. Its per- 
formance improves with higher internal pressures. 

2. The approximate theoretical model developed during this work 
can be used to predict rotating heat pipe performance to within 
+20 percent. 

3- The presence of non-condensible gases can seriously impair the 
heat transfer performance of these devices. 

4. Dropwise condensation gives better performance than film conden- 
sation. This effect is more pronounced at lower pressures and 
rotational speeds. 

5.2 RECOMMENDATIONS 

It is recommended that: 

1. Further work should be performed to improve upon, the heat trans- 
fer characteristics of these devices, including a study of drop- 
wise condensation during centrifugal acceleration. 
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2. In designing a rotating, non-capillary heat pipe, care should be 
taken to ensure that the outside heat transfer coefficient is 
sufficiently high. Otherwise, performance will be impaired. 

3. A study should be performed on the migration and influence of non- 
condensible gases during rotating heat pipe operation. 
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FIGURE 1 SCHEMATIC DRAWING OF ROTATING, NON-CAPILLARY HEAT PIPE 



35 


FIGURE 2 COORDINATE SYSTEM AND GEOMETRY FOR LAMINAR FILM CONDENSATION MODEL, SHOWING 

SCHEMATIC BEHAVIOR OF CONDENSATE FILM UNDER INFLUENCE OF ROTATIONAL ACCELERATION 




FIGURE- 3 
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FIGURE 5 SCHEMATIC DIAGRAM OF CALCULATED FILM THICKNESS 6(x) 
FOR DIFFERENT INITIAL VALUES 5. 
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FIGURE 8 THEORETICAL HEAT TRANSFER RATE VS. SATURATION TEMPERATURE 

FOR WATER SHOWING INFLUENCE OF EXTERNAL HEAT TRANSFER COEFFICIENT 
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HEAT TRANSFER RATE, BTU/HR (X10 


3000 RPM 



FIGURE 9 THEORETICAL HEAT TRANSFER RATE VS. HALF CONE ANGLE FOR 
WATER AT TWO ROTATIONAL SPEEDS 
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FIGURE 10 THEORETICAL HEAT TRANSFER RATE VS. SATURATION TEMPERATURE FOR 
WATER SHOWING INFLUENCE OF INTERNAL CONDENSER RADIUS 
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FIGURE 12. CROSS SECTION OF ROTATING HEAT PIPE 





45 


FIGURE 11. SCHEMATIC DIAGRAM OF EXPERIMENTAL 
EQUIPMENT. 
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FIGURE 13 PHOTOGRAPH OF THE EXPERIMENTAL FACILITY 








FIGURE 15 EVAPORATOR HEAT TRANSFER PERFORMANCE WITH DATA OF OTHER INVESTIGATIONS 
FOR COMPARISON 
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FIGURE 16 CONDENSER HEAT TRANSFER RATE VS. SATURATION TEMPERATURE 
FOR FILM & DROPWISE CONDENSATION 
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FIGURE 17 CONDENSER HEAT TRANSFER RATE VS. SATURATION TEMPERATURE SHOWING 
EFFECT OF FILLING PROCEDURE & PRESENCE OF NON-CONDENSIBLE GASES 
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CONDENSER SURFACE TEMPERATURE 



FIGURE 18 CONDENSER OUTSIDE SURFACE TEMPERATURE PROFILE AT 700 RPM 
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CONDENSER SURFACE TEMPERATURE 



FIGURE 19 CONDENSER OUTSIDE SURFACE TEMPERATURE PROFILE AT 2100 RPM 



FIGURE 20 COMPARISON OF THE CONDENSER EXPERIMENTAL RESULTS WITH THE THEORETICAL 


PREDICTION FOR WATER 
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FIGURE 21 
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COMPARISON OF THE CONDENSER EXPERIMENTAL RESULTS WITH THE THEORETICAL 
PREDICTION FOR WATER 
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FIGURE 22 COMPARISON OF THE CONDENSER EXPERIMENTAL RESULTS WITH THE THEORETICAL 


PREDICTION FOR ETHYL ALCOHOL 
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FIGURE 23 COMPARISON OF THE CONDENSER EXPERIMENTAL RESULTS WITH THE THEORETICAL 
PREDICTION FOR FREON - 113 
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